Finite-temperature phase transitions in a two-dimensional boson Hubbard model 
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We study finite-temperature phase transitions in a two-dimensional boson Hubbard model with 
zero-point quantum fluctuations via Monte Carlo simulations of quantum rotor model, and construct 
the corresponding phase diagram. Compressibility shows a thermally activated gapped behavior in 
the insulating regime. Finite-size scaling of the superfluid stiffness clearly shows the nature of the 
Kosterlitz-Thouless transition. The transition temperature, T c , confirms a scaling relation T c oc p§ 
with x = 1.0. Some evidences of anomalous quantum behavior at low temperatures are presented. 
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Recently quantum phase transitions 1, 2] have drawn a 
lot of attention in systems of interacting particles. Typ- 
ically strong interactions suppress the itineracy of par- 
ticles to induce a strongly correlated insulating phase, 
whereas with weak interactions a conducting phase is 
stable. The criticality of these zero-temperature phase 
transitions can be investigated at low, but finite, tem- 
peratures. How quantum fluctuations associated with a 
quantum critical point(QCP) have influence on phases at 
finite temperatures |3|, |4|, [5| is a theoretically interesting 
and an experimentally relevant question. 

At finite temperatures, it is expected that a quantum 
phase transition turns into a classical one with the same 
order parameter or disappears. Remnant quantum fluc- 
tuations near a QCP may bring anomalous properties [H, 
which can be captured by scaling relations, and lead 
to crossover behaviors as temperature rises. Some pos- 
sibilities such as reentrant behaviors due to the inter- 
play of quantum and thermal fluctuations have been 
proposed[6|. 

These issues can be clarified by direct investigations of 
a generic quantum mechanical model. So far most of the 
theoretical investigations heavily rely on the exact solu- 
tion of the quantum Ising model, available strictly in one 
dimension[4]. Interacting bosonic systems simulated via 
Monte Carlo methods, not suffering from negative sign 
problems, will be an ideal place to study these problems. 
In previous works, a quantum XY model, equivalent to 
hard-core bosons at half-filling, showed the Kosterlitz- 
Thouless(KT) transition^] at finite temperature in two 
dimensions [1, Q. In the model with nearest neighbor 
repulsion, destruction of the solid order as well as the 
superfluidity by thermal fluctuations was observed [loj ] . 
However, generic finite-temperature phase diagrams have 
not been constructed. 

In this work, we investigate the thermally driven phase 
transitions of a two-dimensional quantum rotor model, 
which is believed to share the same critical properties of 
a soft-core generic boson Hubbard model [111], via Monte 
Carlo simulations. The results are summarized in the 
phase diagram as shown in Fig. [TJ Finite-size scaling 
properties of the superfluid stiffness confirm that the na- 



ture of the classical phase transition associated with the 
destruction of superfluidity is consistent with that of the 
KT transition, and clearly support the scenario of the 
universal jump at the critical point 12]. Finite tempera- 
ture, T, sets the size in the temporal direction, leading to 
a scaling behaviorji, 11] T c cx po x with x = 1.0, where T c 
is the transition temperature and po is the superfluid stiff- 
ness at zero temperature. The compressibility diverges at 
the transition. In the insulating regime at low tempera- 
ture, thermally activated behavior of the compressibility 
with a finite energy gap is observed. Some anomalous de- 
pendence of energy and specific heat on T, possibly due 
to quantum fluctuations, are observed for T < 0.251/ . 
The Hamiltonian of a boson Hubbard model reads 



where bj(bA is the boson annihilation(creation) operator 
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FIG. 1: (Color online) Phase diagram on the space of hop- 
ping strength, t(— ^/no(no + and temperature, T, in 
unit of U. The solid line denotes the classical phase transi- 
tions, which terminates at a QCP at T = 0. The dotted line 
represents crossover between gapped fluid and normal fluid. 
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at the j-th site, and rij is the number operator. U and 
w stand for the strengths of the on-site repulsion and of 
the nearest neighbor hopping, respectively, and \x is the 
chemical potential. 

It is convenient to put u/U +1/2 = uq + n with an 
integer no and — 1/2 < n < 1/2 so that no represents the 
background number of bosons per site and fi is a charge 
offset. When n, — 0, the density of bosons is fixed to 
a commensurate filling across the transition. For non- 
integer n, however, an integer filling in a Mott insulator 
shifts to a non-integer filling in a compressible fluid. We 
study the phase transition of the latter case in (2+1)- 
dimensional L x L x L T square lattices, where L denotes 
the size in a spatial dimension and L T in the temporal 
dimension. 

Since the phase transition of the model in Eq. ([T]) 
is characterized by the establishment of phase coher- 
ence, we may rewrite the Hamiltonian in terms of 
the phase angle 6j of bosons by replacing bj(bj) = 
^/nje~^ j (y/rij + le l9j ) with Uj = \-j$p- Under the as- 
sumption that the nature of the transition is governed 
only by the fluctuations of 6j, not those of the hop- 
ping strength, we replace nj — > no so that bj(bj) — 
y/noe~ l9j (y/riQ + le l6j ). Then, the Hamiltonian is re- 
duced to a quantum rotor model 

H = 2 ? ^ ~ 1} ~ 7l] ~ 2t ? ° 0s( - dl ~ 6j) ' {2) 

j j <ij> 

where t = ^Jno{no + l)w. Here we take the number of 
bosons nj > 0. 

Through a path integral map ping , we can construct 
the corresponding classical action [14j 

S = ^2 ^- Jr(Jr ~ 1) - e M- - In Jjj. (2et) - Inljy (2et)(3) 

r 

with the partition function 

V-,7=0 

Z= E e~^\ (4) 
{f r } 

where e = /3/L T is a lattice constant in the imaginary 
time axis for an inverse temperature 0, J r is an integer 
current at site r = [j, r) with a spatial index j and a tem- 
poral index r, which is conserved at each site as denoted 
by V • J = 0, and I m (x) is the modified Bessel func- 
tion given by the relation e Kcos9 = YZ=-oo I m {K)e lm6 . 
In this work, we investigate the properties of the model 
in Eq. [3] via Monte Carlo simulations using a recently 
proposed worm algorithm [l3| . In order to reduce the 
systematic errors in discretizing the imaginary time axis, 
we need to take eVtU -C 1. We take Ue = 0.5-2 for 
t C 17 and set the energy unit {7 = 1. 

The superfluid stiffness in a finite system is given bvflij] 

p L =p-*I?-*(WZ), (5) 
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FIG. 2: (Color online) Finite-size scaling behaviors of the su- 
perfluid stiffness as a function of (a) hopping strength and 
(b) temperature. For both cases, data collapsing onto a sin- 
gle curve works fine in terms of the scaling parameter L/£ as 
shown in insets, consistent with the nature of the KT transi- 
tion and the universal jump at the critical point. 

where W x = L^ 1 ^ jr J^ and (...) denotes the averages 
over the probabilites determined by the partition func- 
tion of Eq. ([¥]), and d is the spatial dimensionality. Sim- 
ilarly the compressibility is 

K = f3L- d {(N 2 ) - (N)% (6) 

with N = L~ 1 ^ jr . The energy expectation is given 
by 

(H)=L-\—), (7) 
oe 

and the specific heat is C v = L~ d (d{H) / dT). 

We consider the case for /j, = 0.9 so that no — 1 and 
n = 0.4. Figure [2] shows the finite-size scaling behav- 
ior of the superfluid stiffness as a function of (a) t and 
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(b) (3. Finite-size scaling properties of the transition can 
be obtained by plotting the curves in terms of a scal- 
ing variable £/£, where £ is the correlation length. Here 
we assume an essential sing ularitypjl £ ~ exp^S' 1 / 2 ), 
where 8 = t — t c (or (3 — f3 c ) is a tuning parameter and 
& is a non-universal scaling factor. In terms of this scal- 
ing variable, we obtain high-quality data collapsing onto 
a single curve for different sizes, consistent with the na- 
ture of the KT transition. The scaling behavior also sup- 
ports the scenario of the universal jump of the superfluid 
stiffness 12J, (7r/2)/3 c p 00 = 1, at the critical point in the 
thermodynamic limit. By extrapolating the single curves 
to the critical point, we find that (w/2)/3 c p 00 w (a)1.01 
and (b)1.06. These numbers are, however, sensitive to 
fitting parameters b and t c ((3 c ). 

Figure [3k shows the behavior of the compressibility. 
The finite-size scaling ansatz of the compressibility is 
written in the form 



K = L z - d X K {L{t-t c ) 1 !",p/L*), 



(8) 



where X K is a dimensionless scaling function and z is the 
dynamical critical exponent. For the generic superfluid- 
insulator transition(GSIT), z = 2 is expected [111]. The 
crossing behavior of the compressibility curves for differ- 
ent sizes at the critical point t° c = 0.023±0.001, therefore, 
represents the scaling properties near the QCP, where i° 
is the critical hopping strength at zero temperature. For 
different values of p, we have similar results with tj. just 
shifted. 

We find that the compressibility diverges at the transi- 
tion. In the superfluid side, ft ~ l/(t — This strongly 
supports that the longe-range density fluctuations drive 
the transition. In the insulating side, the compressibility 
has an activated form e~ Agap / T with a finite energy gap 
Agap- This dependence is shown in Fig. [3)3 for different 
t, from which we can calculate A gap as shown in the in- 
set. For small t, we need a large number of Monte Carlo 
steps to obtain equilibrium and have bigger error bars in 
determination of A gap . The gap vanishes around t = t® 
as expected. 

Thus we have a so-called 'V-shaped' phase diagram 
(FigQ]). In the insulating regime, the Mott insulator 
exists at T = 0, which turns into an activated gapped 
fluid with a finite energy gap at low temperature. It 
gradually disappears in a high-temperature normal fluid. 
This crossover line can be specified by the condition 
Agap/T w 1. The phase coherence in a superfluid at 
T = is destroyed by quantum fluctuations to form a 
QCP or by thermal fluctuations at T > to define clas- 
sical phase transitions. The phase boundary in Fig. Q] is 
obtained by tuning t for given T (black circles) as well 
as by tuning f3 for a given t (red squares) . Note that the 
phase boundary follows a scaling relation T c oc \t — t® c \ zu , 
which implies that /3 determines the correlation length in 
the temporal direction, where v is the correlation length 



L=12 L =18 

T 

L=16 L t =32 
L=20 L t =50 - 
L=24 L t =72 
L=28 L =98 





1=005 



P 



FIG. 3: (Color online) (a) Compressibility of the boson Hub- 
bard model shows behavior of the GSIT with z — 2.0, diverg- 
ing at the transition, (b) In the insulating regime, we have 
thermally activated behaviors, k ~ e - A s*p/ T j from which 



Ag a p can be evaluated, 
ishing at the QCP. 



Inset: 



critical exponent. The boundary in Fig. [T] is consistent 
with the expectation zv = 1 [1 1] for the GSIT. 

It is interesting to check the predicted scaling relation 
0, EH T c oc p x in this model. Figure S] shows that 
the zero-temperature superfluid stiffness p°, denoted by 
dotted line, which obtained via extrapolation of values at 
T > 0, follows p° cx \t - t% implying that x = 1.0. It is 
consistent with the hyperscaling argument 11] suggesting 
x = z/(d+ z - 2). 

We expect that this quantum criticality disappears 
as temperature rises, which means quantum fluctuations 
possibly leave some tracks in bulk properties at low tem- 
peratures. Figure [5] shows the specific heat, Cy, and the 
energy expectation values, (H), as a function of T for dif- 
ferent t. Sharp rises of Cy in the conducting regime or 
round up-rises in the insulating regime are followed by 
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indents, regions indicated by A, which apparently rep- 
resent anomalous behavior due to quantum fluctuations 
and disappear at high temperatures for T > 0.25. This 
feature strongly suggests a crossover in normal fluid from 
quantum mechanical to classical regime. Similarly the 
curves of (H) show bumps, indicated by T, only in the 
range where quantum critical fluctuations are expected 
to have effects. 

In summary, we have investigated the phase transi- 
tions at finite temperature in a two-dimensional quan- 
tum rotor model in which intrinsic zero-point fluctua- 
tions are present. Finite-size scaling of the superfluid 
stiffness shows an essential singularity of the KT phase 
transition and the universal jump at the critical point. 
The compressibility diverges at the transition. In the 
insulating regime, the compressibility shows a thermally 
activated behavior, k ~ e - A ga P /^ f rom wn ich we can 
successfully evaluate the gap. This indicates that the in- 
sulating behavior at low temperature gradually crosses 
over to the behavior of normal fluid as temperature in- 
creases. The transition temperature T c shows a scaling 
behavior T c oc \t — t°|, showing that finite T limits the 
length of quantum fluctuations in the temporal direction, 
and a hyperscaling relation T c oc po- The behavior of the 
specific heat and the energy suggests that, as tempera- 
ture rises, quantum critical regime near a QCP crosses 
over to classical regime. 
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FIG. 4: (Color online) Superfluid stiffness for different j3. As 
13 increases, the size dependence becomes smaller. This allows 
us to extrapolate the curves to obtain zero-temperature su- 
perfluid stiffness, po, in the thermodynamic limit as denoted 
by dotted line. It shows that po oc \t — t°\ with w 0.22. 
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FIG. 5: (Color online) Specific heat, C'v, as a function of T 
for different t. Sharp rises in the conducting regime, signa- 
ture of the superfluid transition, or round up-rise of Cv in the 
insulating regime are followed by indents which disappear in 
high temperature region, T > 0.25. Insets: The curves of the 
energy expectation values, {H), have bumps at low tempera- 
tures possibly due to the effects of quantum fluctuations. 



